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A model for deeply virtual Compton scattering, based on analytical light-cone 
hadron wave functions is presented and studied at energies currently accessible at 
Jefferson Laboratory and DESY. It is shown that poles and perpendicular vector 
components play an important role at < 10 GeV'^. A suppressed diagram 
has to be included at these low energies, but becomes negligible above 10 GeV^. 
Future prospects and developments of this model are discussed. 



1. Introduction 

For many years, deeply inelastic scattering (DIS) has been a major source 
of our knowledge of parton distributions in nucleons and nuclei. This is 
because it can be shown that, at sufficiently high energies and momentum 
transfers, the amplitude for this process factorizes into a 'hard' part which 
can be calculated from QCD, and a 'soft' part which can be extracted 
from experimental data. The soft part can be proved to be related to 
the probability of finding a quark with a particular flavor carrying a given 
fraction of the nucleon's momentum. 

In recent years, much interest has been focused on studying deeply vir- 
tual Compton scattering and the electroproduction of mesons, following 
the proof by Collins, Frankfurt and Strikman^. This proof demonstrated 
that, under quite general conditions, the leading amplitude for hard exclu- 
sive photo-production of mesons could also be factorized into a calculable 
'hard' part, and a 'soft' part. All other amplitudes were smaller than the 
leading amplitude by powers of l/Q. The soft part corresponds to the pro- 
cess by which a parton with a certain momentum fraction is removed from 
a nucleon, and replaced by a parton with a different momentum fraction. 
These parton distributions have been given a variety of names, but we will 
refer to them as 'generalized parton distributions' (GPD's). In the case of 
quark helicity conservation, there are four independent GPD's; H{x,(,t), 
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E{xX,t), H{x,(,t), and E{x,(,t), where x and ( are the hght-cone mo- 
mentum fractions of the struck quark and real photon, while t = is the 
momentum transfer squared. Three physically different regions could be 
distinguished for x and (. The domain 0<C<2;<1(C^1<2;<0) 
corresponds to the removal and replacement of a quark (antiquark) with 
momentum fractions x {( — x) and x — ( {—x), respectively. In the remain- 
ing region < a; < ^, the photon scatters on a virtual quark-antiquark pair, 
extracted from the proton. With this notation, C ^ xb (Bjorkcn x) in the 
limit ^ oo {A fixed). In the limit of forward scattering (DIS), the ff's 
reduce to the quark density and quark helicity distributions: 

H{x,0,0) = q{x) (1) 
H{x,0,0) = Aq{x). (2) 

The E^s do not appear in DIS, they are unique to the off-forward exclusive 

processes and provide information not accessible through other means. The 
GPD's are related to the nucleon form factors by the integrals 

dx 
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^_^_m{x,C,t)=F^{t), (3) 
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E^{x,C,t)=F^{t), (4) 
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and similarly for H and E, all independent of C- The factor 1 — | is included 

to comply with the normalization of form factors used by Ji^. The GPD's 
might also shed some light on the nucleon spin decomposition, since Ji'^ has 
related them to the quark spin: 

1 dr 

-^"=0 —^x[m{x,C,t = 0)+E'i{x,C,t = 0)]. (5) 

^ ^c-i 1 - I 

The quark angular momentum is decomposed into Jq = Sq+Lq, where Sq is 
measured in polarized DIS experiments. If the above sum-rule is measured 
in DVCS, the quark spin and orbital angular momentum parts could be 
separated. 

An alternative approach to DVCS is offered by the light-cone quark 
wave functions as suggested by Brodsky, Diehl, and Hwang"^ (BDH). They 
have used this idea to calculate DVCS on an electron in QED for large Q^. 

There are ambitious efforts under way to measure DVCS (and the hard 
exclusive meson photo-production) at a number of facilities, DESY and 
JLab in particular'*. At JLab energies, the competing Bremsstrahlung or 
Bethe-Heitler (BH) process is larger than DVCS. However, by doing in- 
terference measurements (e^/e" beam charge asymmetry and various spin 



Figure 1. The five diflferent time-ordered 'handbag' diagrams. 



asymmetries) the BH amplitude cancels out and only a BHxDVCS inter- 
ference remains. At DESY the energy is large enough for DVCS to become 
larger than BH, though the measurements have low statistics, not allow- 
ing for differential cross sections. Data for the asymmetries are, however, 
at both facilities collected for relative low beam energies and momentum 
transfers. It is thus interesting to investigate to what extent the 'leading' 
amplitude is actually dominating in this kinematic region. We will here 
present the first few steps toward such an understanding, by developing a 
simple model using effective analytic quark-diquark wave functions, which 
allows us to study the amplitudes excited at these low energies. This model 
is similar to the one of BDH, but we calculate DVCS on a proton and focus 
especially on the features at low Q^, keeping higher twist terms. 

2. Formalism 

We start from the Fourier transform of the 7*p 7p' amplitude (in light- 
cone coordinates) 



where J^{y) = '4'{y)l^ Ib^i^iv) is the electromagnetic current and p{p') 
and q{q') are the four-momenta of the initial (final) hadron and photon. 
This expression could be expanded to give the five different one-loop dia- 
grams (assuming light-cone gauge = 0) shown in Fig. 1. These diagrams 
are essentially the same as those used by BDH for DVCS on an electron. 
The diagrams are represented by the integrals (because of space limitations, 
only a and b are given) 
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where z = {x — — 0, z' = {(^ — x)/{l — x), and the relative momenta 

are defined as 1_l = k_L + (1 — z)A±, 1^ = —(1 — z')k± — A±, = 
k_L + (1 -x)q_L, and in±b = k_L + + ^^q_L. Here x = k'^/p'^ > and 
C = q'~^ /p'^ are the longitudinal momentum fractions of the struck quark 
and the real photon, m, mji, and M the quark, remnant (diquark), and 
hadron masses, and A = q — q' ■ The expressions for the other diagrams 
have similar structures. In the limit — > oo, the denominator of Eq. (8) is 
proportional to x — i^ + ie, i.e., the leading twist expression of Ji^. Note that 
the integrals are over different parts of phase space and include scattering 
on a quark or quark-antiquark pair only, not scattering on an antiquark. 
The wave function is represented by the analytical form 



{x, kj_) = iVexp 
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where (3 = 0.69 GeV is chosen such that F{q^ 
kj_)(/)(x,kj_) agrees with the dipole form factor for < 
skew diagrams require the knowledge of the wave function (p{z',l'j_) for 
the diquark splitting into a hadron and quark (lower right-hand corner of 
diagrams b, c, and e). The form of this wave function will eventually be 
restricted by exclusive data, but is here arbitrarily chosen to be of the form 
of Eq. (9), with mR and x replaced by M and z'. 



3. Results 

The r++ matrix elements have been calculated for present JLab and DESY 
kinematics and are shown in Fig. 2. The laboratory angle 9^^' between the 
virtual and real photons is defined for in-plane kinematics such that it is 
positive for cf) = and negative for cj) = 180°, where cj) is the azimiith angle 
between the final electron and proton, with q as the polar axis. These 
calculations incorporates the principal value parts of the integrals, i.e., the 
imaginary part of the T++. 

The forward peaking of the skew diagrams (b, c, and e) is a consequence 
of the shift of the momentum kj^ in the extra denominator. The maximum 
of the integral is dislocated to momenta where this denominator becomes 
very small. 
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Figure 2. The ImT++ matrix elements as functions of the laboratory angle between 
the virtual and the real photon for in-plane scattering. The labels correspond to the 
ones of Fig. 1 (the lower solid line always correspond to diagram d). The calculations 
are for xb = 0.35 and = 2,4, 10,20 GeV^ from top to bottom. 

Because of its two hard propagator, diagram b is expected to be sup- 
pressed in the high hmit^, and our calculation indicates that this sup- 
pression becomes significant for Q"^ > 10 GeV^, i.e., well above present 
JLab energies. 

While diagrams a and b exhibit a physical pole (for massless final pho- 
ton), the other diagrams never get vanishing denominators for physical val- 
ues of kinematic parameters. Thus the real part of the DVCS amplitude, 
which comes from the i5-function part of the propagators, has contributions 
from these first two diagrams only. The real part calculations will not be 
presented here. 
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This difference in pole structure also explains why diagram d is smaller 
than diagram a, since terms add up in the denominator of d, while they 
could cancel for a. This feature is closely related to the behavior of the two 
leading-twist propagators of Ji^; l/{x — ^ + ie) and + ^ — ie), where x 
and ^ are momentum fractions related to ^{p~^ +p'^) instead of p+. In this 
notation there is a cancellation for a; = ^ (scattering on quark) or a; = — ^ 
(antiquark) . These simplified propagators neglect four- vector components, 
e.g., A_L, that do not give large scalars in the Bjorken limit. 

4. Conclusions cind outlook 

We have introduced a model for DVCS, using simple quark wave functions 
and retaining all components of four-vectors throughout. The perpendicu- 
lar momenta are shown to be very important for the location of the maxima 
of the integrals at low Q^. Within the model we have been able to investi- 
gate the five single-loop diagrams and their relative importance for DVCS 
at JLab and HERA energies. The results indicate that at present JLab 
energies (Q^ < 4 GeV^), all of the diagrams (except possibly the crossed 
diagram d) need to be considered. In particular it is necessary to include 
diagram b, despite its two hard propagators, since the expected suppres- 
sion is significant only for larger Q^. Other higher-order diagrams might 
need to be considered as well. At > 10 GeV^, the process is completely 
dominated by the handbag diagrams with one hard propagator. 

This work will be extended to cover a wider kinematic range in 
and Xb and will include antiquark contributions to the amplitudes and the 
wave functions. The photo-production of mesons will also be calculated, 
with special consideration of the meson poles then made possible in the 
skew diagrams. 

In order to compare with actual JLab and HERA experiments the inter- 
ference with the Bremsstrahlung (Bethe-Heitlcr) process has to be calcu- 
lated. We would then obtain cross sections and asymmetries to test against 
experimental results and be able to check the validity of commonly used 
approximations and to evaluate various wave functions that could be used. 
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